Scale-invariance of human EEG signals in sleep 
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We investigate the dynamical properties of electroencephalogram (EEG) signals of human in sleep. 
By using a modified random walk method, We demonstrate that the scale-invariance is embedded 
in EEG signals after a detrending procedure. Further more, we study the dynamical evolution 
of probability density function (PDF) of the detrended EEG signals by nonextensive statistical 
modeling. It displays scale-independent property, which is markedly different from the turbulent- 
like scale-dependent PDF evolution. 

PACS numbers: 87.19.Nn, 05.40.-a, 87.20.-a, 89.75.Da 



-Introduction. The analysis of electroencephalogram 
(EEG) signals attracts extensive attentions from various 
research fields, since it can not only help us to under- 
stand the dynamical mechanism of human brain activi- 
ties, but also be potentially useful in clinics as a criterion 
of some neural diseases. Some previous works have been 
done on human EEG signals in sleep and other physio- 
logical states. In Refs. [J 12, 13] the correlation dimension 
and Lyapunov exponent are calculated to characterize 
and discriminate the sleep stage. Lee et al. provides the 
evidence of the long-range power law correlations embed- 
ded in EEG signals The mean scaling exponents are 
distinguished according to REM (Rapid Eye Movement), 
Non-REM and awake stage, and gradually increased from 
stage 1 to stage 2, 3 and 4 in non-REM sleep. Hwa et 
al. found the variable scaling behavior in two regions, 
and described the topology plots of scaling exponents in 
this two regions that reveals the spatial structure of the 
nonlinear electric activity The random matrix the- 
ory is performed to demonstrate the existence of generic 
and subject-independent features of the ensemble of cor- 
relation matrix extracted from human EEG signals Q. 
Yuan et al. found the similar long-range temporal cor- 
relations and power-law distribution of the increment of 
EEG signals after filtering out the a and /3 wave @- In 
the present paper, the Tsallis entropy is used to analyze 
a series of human EEG signals in sleep. 
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We use the MIT-BIH polysomnography data, which is 
consist of four-, six- and seven-channel polysomnographic 
recordings, each with an ECG signal annotated beat-by- 
beat, and an EEG signal annotated with respect to sleep 
stages Q. Records have been sampled at frequency 4 
kHZ. Sleep stage was annotated at 30s intervals according 
to the criteria of Rechtschaffen and Kales, denoted by six 
discrete levels-1, 2, 3, 4 REM and awake (stages 1, 2, 3, 
4 belong to non-REM sleep) [9]. In the present analysis, 
only the samples containing sufficient records (at least no 
least than five stages) are considered. A representative 
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FIG. 1: A set of representative records of EEG signals in 
different stages. Each entire series includes more than 10 5 
data points, while the plotted is only a small fraction. 



example is shown in Fig.l. 

-Scale-invariance of detrended EEG signals. Consider 
an EEG series, denoted by {xi] (i = 1, 2, . . . , N), whose 
scaling characteristics are detected through the following 
procedure: 

Step 1: Construct the profile series, Yj — Yll—i Xi,j = 
1, 2, • ■ • , N, and consider Yj as the "walk dsiplacement" 
of a resultant random walk. 

Step 2: Divide the profile series into non-overlapping 
segments with equal length and fitting each segment with 
a second order polynomial function. Regard the fitting 
results as the trends, a stationary series can be obtained 
by eliminating the trends from the profile series. 

Step 3: After the detrending procedure, we define the 
increment of this modified profile series at a scale s as 
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where Y* is the deviation from the 
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FIG. 2: The probability of P(0, s) as a function of the time 
sampling scale s. A power-law scaling behavior is observed 
for about three order of magnitude. 



FIG. 3: (color online) The rescaled increment's PDFs for six 
stages. Obviously, curves with different time scales can well 
collapse into a single master one, demonstrating the existence 
of quasi-stationary property. 



polynomial fit. 

step 4 : Scale-invariance (self-similarity) in the sta- 
tionary series implies that the PDF satisfies, P(x, s) = 
■£-P {.■§-) i where a s denotes the standard deviation at 
time scale s. Obviously, P(0, s) — P(0) — . 

Assigning the values of parameter s from 2 1 to 2 10 , 
the normalized PDFs of A S Y exhibit scale-invariant (self- 
similar) behaviors as presented in Fig. 2. That is to say, 
those PDFs can be rescaled into a single master curve, as 
shown in Fig. 3. The scale-invariance of the detrended 
EEG signals suggests that the quasi-stationary property 
is embedded in the distributions of time scales. There- 
fore, it helps us to search for stable distributions to mimic 
them. 

-Nonextensive statistical modeling of detrended EEG 
signals. From the results sketched in the preceding sec- 
tion, herein we use the Tsallis entropy to model the 
PDFs. The Tsallis entropy is introduced by Tsallis 
through generalizing the standard Boltzmann-Gibbs the- 
ory [1(1, which reads, 
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1 Jdx\p{x)]« { [ dxp{x) = 1;qeR) . (1) 
9-1 J 

In the limit q — ► 1, S q degenerates to the 
usual Boltzmann-Gibbs-Shannon entropy as Si — 
— J p(x)ln[p(x)]dx. The optimization (e. g. maximize 
S g if q > 0, and minimize S q if q < 0) of S q with the nor- 
malization condition J dxp(x) = 1, as well as with the 



FIG. 4: (color online) The rescaled increment's PDFs of all 
stages with the approximate fit using nonextensive statistical 
modeling. We use g-Gaussian function to fit the awake stage, 
and g-exponential function to fit the other five stages. The 
awake stage falls into levy regime with well fitting parameter 
q — 1.94. In the REM stage, the values of q is slight change 
describing as q2, ■ ■ ■ , 91024, while in the non-REM stages, 
parameter q shows a unique value for each specific stages. 
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value. 

The results in Fig. 4 show us that the PDFs of awake 
stage falls into the Levy regime with q being equal to 
1.94. It exhibits sharp kurtosis and long-tail distribution, 
distinguished from those of REM and non-REM stages. 
It should be noted that we shift the distributions through 
dividing by their standard deviation and only plot part 
of them to make the figures clear. The specific values of 
/3(s) for all scales are shown in the Fig. 5. It interests us 
that /3(s) does not dissipate as the scale increases unlike 
the case of (3(s) found in financial market In other 
words, It demonstrates that the dynamics evolution of 
EEG signals is not coincident with the diffusion process 
described by Fokker-Planck equation. 

Another significant equation of nonextensive statistical 
approach is the g-exponential function, which reads 



FIG. 5: (color online) /3(s) and r(s) versus s for awake and 
non-REM stages. The values of /3(s) don't dissipate as the 
increasing of s. In particular, r(s) of non-REM sleep converge 
to an invariant pattern. 
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FIG. 6: (color online) The increment's PDF of randomized se- 
ries of awake stage and fitting curves with different parameter 
q. The parameter q rapidly approaches to Gaussian regime 
(q = 1) as the time scale increases. For clarity, we shifted the 
distributions through dividing them by their standard devia- 
tion 

constrain ((x 2 )) q = er 2 , leads to g-Gaussian distribution 

= ^t{1-^(s)[(1-9)(^-S(s)) 2 )]}^, (q < 3) 

(2) 

where Z q (s) is a normalization constant, /3(s) is explic- 
itly related to the variance of distribution, and the sub- 
script " + " indicates that G q (x, s) is non-negative [Til ]. 
G q ->i(x, s) recovers the usual Gaussian distribution. The 
q-Gaussian PDF can describe a set of stable distribu- 
tions from Gaussian to Levy regimes [l2| by adjusting 
the value of q with appropriate time-dependent parame- 
ters 0(s) and Z q {s) [lj]. The distribution falls into Levy 
regimes in the interval | < q < 3, with q — 5/3 a critical 



e q (x, s) = ——{I - t(s) [(1 -q)\x- x(s)\}}^ , (3) 
^ q (s) 

where the parameter r(s) is the relaxation rate of dis- 
tribution. Clearly, in the limit q — > 1, e\{x,s) = 
z 1 ^ exp{~r(s)\x — x(s)\}. Because the statistic distri- 
butions of detrended increment of EEG signals in sleep 
stage exhibit an approximately exponential form, we use 
the g-exponential model to quantifies them, as shown in 
Fig. 4. The values of q for the REM and non-REM stage 
are little larger than 1. It means that the fluctuation of 
human brain activity in sleep stage will converge to a nor- 
mal exponential pattern. In particular, the EEG signals 
exhibit g-invariant pattern for different time scales in all 
the four stages within non-REM sleep. The relaxation 
rates of distributions are also approximately invariant, 
as shown in the Fig. 5. However, in the REM stage, the 
values of q are slight change because of fitting the tail 
of distribution in different time scales, and the model 
can only well fit the center distribution. It suggests that 
brain electric activity in the REM stage may work in a 
more complex pattern than awake and non-REM stage 
for the acute neural activity (l5| . 

The nonextensive statistical approach modeling the de- 
trended increment's PDF of EEG signals with an in- 
variant parameter q demonstrates the scale-independent 
property of the system. In order to further test the exis- 
tence of this observed property, we randomize the empir- 
ical series of awake stage by shuffling [I?} and show a 
fit for this artificial distributions at different scales in the 
' Fig. 6. Clearly, the parameter q will approach to Gaus- 
sian regime (q = 1) as the time scale increases. This 
result strongly supports the scale-independent property 
of human brain activity in sleep is remarkably different 
from the turbulent-like scale-dependent evolution [lSj |. 

-Conclusion. In this work, several dynamical proper- 
ties of human EEG signals in sleep are investigated. We 
firstly use a modified random walk method to construct 
the profile series including the information of EEG sig- 
nals. After a detrending procedure, we obtain the sta- 
tionary series and define the increments of the resultant 
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random walk at multiple scales. In order to character- 
ize the dynamical process of brain electronic activity, we 
then study the P(0, s) of the PDF of normalized incre- 
ments as a function of s. With this choice we investigate 
the point of each probability distribution that is least ef- 
fected by the noise introduced by experimental data set. 
The scale-invariance in both awake and sleep stages are 
obtained, thus one can rescale the distributions at differ- 
ent scales into a single master curve. 

Aim at this property, we use nonextensive statistical 
approach to model these processes. The dynamical evo- 
lution of dctrendcd increment's PDF in awake stage can 
be well fitted by g-Gaussian distribution with a invariant 
parameter q — 1.94. It demonstrates that the PDFs of 
awake stage fall into the levy regime. Contrastively, q- 
exponential distribution is used to mimic PDFs of sleep 
stage. In particular, for the non-REM stage, it exhibits 
scale-independent distributions, while for REM stage, it 
suggests a complex distributional form with the values of 
q slightly different. 

The statistical properties of distribution of EEG sig- 
nals strongly indicate that the process of brain electric 
activity is remarkably different from tubulent-like cas- 
cade evolution. In a recent work [18(, Lin-Hughson pro- 
posed a tubulent-like cascade model to mimic the human 
heart-rate, whose validity is, now, in the face of challenge 
on critical scaling- invariance found in real human heart- 
rate processes [ill Ho[ ■ Although the electrocardiograph 



(ECG) and electroencephalogram are different, they have 
some common statistical features. We wish this work is 
helpful for the in-depth understanding about the under- 
lying dynamical mechanism. And, like the correspond- 
ing empirical studies on human ECG singles [19|,l20(, this 
work could provide some criterions for theoretical models 
on human EEG signals. 

These similar results are also found in the ECG signals 
of human. Since the work presents the turbulent-like cas- 
cade heart rate model to mimic human heart rate, Kiyono 
et al. found that the human heart rate exhibits criti- 
cal scaling-invariance, of which the dynamical evolution 
of increment's PDF is different from turbulent-like PDF 
evolution. 

Further more, some very recent works [22I |23| pointed 
out that the sleep-wake transitions exhibit a scale- 
invariant patterns and embed a self-organized criticality 
(see also Ref. [2l[ for the concept of self-organized crit- 
icality). Thus, the dynamical properties of human EEG 
signals in sleep suggest that human brain activit y in sleep 
may relate with a self-organized critical system [24| . Our 
empirical result, in some extent, support those conclu- 
sions. 

This work is supported by the National Natural Sci- 
ence Foundation of China under Grant Nos. 70571075, 
70571074 and 10635040, and the Foundation for graduate 
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